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ABSTRACT 

We present the result of IR observations in the K band of 39 bright AeBe stars in the 
northern hemisphere. Using the shift-and-add technique, we were able to detect nearby objects 
as close as 0.4" from the primary star. We found 9 stars that have one or more companions. 
Adjusting for completeness, we compute that the AeBe star binary rate is 85%. The binary rate 
of AeBe stars is therefore greater than that of near solar type MS stars (57%) and similar to 
that found in T Tauri stars (80%). 

Subject headings: infrared:stars — stars:pre-main-sequence — stars:statistics 



INTRODUCTION 



AeBe stars are intermediate mass pre-main sequence (PMS) stars (1.5 < M /Mq < 10) first recognized by 



Herbig (1960) to be the higher mass equivalents to the T Tauri stars (TTS). Herbig (1960) defined AeBe 



stars to be stars of spectral type A or earlier, located in an obscured region, and illuminating a bright 
nebulosity. The number of stars in Herbig's original list was 26, a number later expanded to 57 stars by 



Finkcnzcller fc Mundt (1984) . While originally broad, the definition of an AeBe star has been steadily 



refined over the years (see The et al. (1994) for a good discussion) as the number of AeBe candidates grew. 
To date, the most comprehensive list of AeBe stars can be found in the work of The et al. (1994)| who 
published a list of more than 287 AeBe candidates. 



After an extensive survey of 47 stars Hillenbrand et al. (1994) distinguished among three distinct groups 
of AeBe stars: stars with large infrared excesses (Group I, candidate disk objects), stars with a flat or 
rising spectra (Group II, envelope/disk objects), and stars with modest infrared emission (Group III, 



diskless objects). Hillenbrand et al. (1994) proposed that the infrared excesses seen in Group I AeBe stars 
were the result of circumstellar disks. They showed that the nature of the infrared excesses (specifically a 
spectral energy distribution (SED) of the form \Fx = A^3 ) could be well modeled by a centrally evacuated 
(3 < Rhoie/R* < 25) optically thick circumstellar disk. 



Hartmann et al. (1993) later argued that the high accretion rates required in the Hillenbrand et al. 
model to account for the observed amount of infrared light were incompatible with the existence of an 
evacuated region at the center of the disk. Instead, relatively small IR excesses could be generated by small 
dust grains heated in a reflection nebula surrounding the central star, while larger IR excesses might be 



produced by close stellar companions. The importance of close nearby objects was explored by Kenyon et 



al. (1993)1 , who successfully modeled 1-3/im infrared light excesses as being produced by deeply embedded 
stars. 



To help to resolve this issue, Li et al. (1994) attempted to determine the binarity of AeBe stars by 
obtaining quasi-simultaneous JHK photometry of 19 AeBe stars. They used small apertures to avoid 
contamination from nearby sources and found that 9 stars had a companion star within 10 arcscconds 
of the primary with separations ranging from 2500 to 9000 AU. Their search was however restricted by 
atmospheric seeing to objects farther than 1.1" - 1.6". For an AeBe star at a distance of 2500 pc, such as 



MWC 1080, sources closer than 2750 AU to 4000 AU would not have been detectable. 



High resolution observations can be achieved from a ground based telescope using the shift-and-add 



(SAA) technique (Bates and Cady (1980)). SAA is a simple technique which consists of taking a series of 
very short exposures (specklegrams) that freeze the image of the star before it is blurred across the detector. 
A long exposure image, free of atmospheric tilt errors, can then be constructed by individually shifting the 
image of the star to the middle of each specklegram before co-adding them together. This technique has 
been shown to work well in the past and to retrieve most of the high spatial frequencies usually lost to 



seeing ( [Christou (1991)| ) 



Our experience with SAA at a 2.3m class telescope has shown that we are typically able to reduce the 
FWHM of a point source by a factor of 3 in the near infrared. On observing nights with an average seeing 
we repeatedly managed to construct high resolution images with a FWHM of 0.3" - 0.4". 

We conservatively estimate our inner detection limit at 0.4" and maximum K magnitude for an observable 
secondary at K = 10.5 (with S/N of 10 in the final shift-and-added frame). This is to be compared to the 
angular resolution limit by Rayleigh's criterion of 0.24". 

2. OBSERVATION 

The data presented here were obtained in July and December 1995 and July 1996 using the WIRO 2.3m 
telescope and a 256x256 NICM0S3 infrared camera equipped with the approximate Johnson J, H, and K 
filters. The image scale was .12" /pixel and we only used one quadrant of the detector at a time resulting in 
a field of view of 15". In July and December 1995, each observation consisted of a series of one hundred 60 
ms specklegrams obtained in less than a minute, thus ensuring that the point spread fmiction (PSF) of the 
system did not change during the course of the observing. The data were taken in the same manner in July 
1996 but using larger series of specklegrams, each containing 200 frames. 

Due to our short integration times, the readout noise of the camera was the dominant source of noise. 
Hence, since we required a S/N ratio high enough to accurately determine the position of the primary star 
in each specklegram, we were only able to observe objects brighter than 9*'* magnitude (in the K band). 



The S/N of secondary objects in the field of view did not have to be as high as demonstrated by the fact 
that we were able to detect companions as faint as 10.5 magnitude. 



We observed 39 of the brightest AeBe stars listed by The et al. (1994) over a period of one year in the K 
band (Table 1). Each program star was preceded and immediately followed by a point source star to serve 
as an estimate of the PSF. This was done twice using different quadrants of the detector. The specklegram 
series was then shift-and-added to produce high resolution images of the program stars and point source 
objects. 

3. DISCUSSION 

We were able to identify 9 multiple systems (7 binaries and 2 triple systems) with angular separations 
ranging from .6" to 7" (figure 1). For statistical reasons we differentiate between two types of stars: those 
with a reasonable distance estimate (27 stars), and those with no distance estimate (12 stars). Among the 
first group we found 7 multiple systems with projected physical separations ranging from 176 to 11725 AU. 
Two additional multiple systems were found in the second group (HD 216629 and HD 141569). Previously 



observed companions around HD 200775 (Li et al. (1994)), and around XY Per (Herbig & Robin-Bel] 



(1988)) were successfully detected. The separation and position angle of the detected multiple systems are 



given in table 2. The close companions around MWC 1080, XY Per, and HD 150193 were successfully 



resolved using the CLEAN algorithm ( Hogbom (1974) ), using the image of a nearby star to serve as an 



estimate of the PSF. The two companions surrounding XY Per are an artifact from the shift-and-add 
caused by the fact that XY Per and its companion are close and are of nearly equal brightness. 

Using these observations, it is possible to determine the mean number of companions for AeBe stars and 



to compare it to that of TTS (Ghcz et al. (1993)) and near solar type MS stars (Duquennoy & Mayor 
(1991)). For each object in our list, a Monte Carlo technique was used to generate and project onto the 
background sky ten thousand binary systems with all possible orientations (i, w, 17) and configurations (a, 
e, q—M2/Mi). These were used to determine the number of times a companion would be visible around a 
given AeBe star. In order to be visible, a nearby companion had to be between 0.4" and 8" (conservative 
estimates of the resolution limit and of half of the field of view, respectively) from the primary star and 



had to be massive enough to be visible in the K band (we assumed main sequence color corrections, Allen 



(1973)). The eccentricity and mass ratio distributions were taken from Duqucnnoy & Mayor (1991). The 
masses of the AeBe primaries were not known and we averaged the results obtained using a flat distribution 
of masses ranging from l.SM© to IOMq. When the distance to an AeBe star was unknown, we additionally 
averaged the results obtained by placing the star at distances ranging from lOOpc and ISOOpc. Table 1 
shows the thus predicted binary discovery rate of each of the stars in our sample. 

By assuming an intrinsic binary frequency of 100%, we computed that 16% ±1.6% of the companions 
would be detected or 6.2 ±0.6 stars for a sample of 39 binary stars. We also calculated that 0.2 field stars 



brighter than 10.5 K magnitude would be observed within 8" (based on computations using the Mamon 



fc Soneira (1982)| model). Adjusting for fortuitous field stars, we therefore expected to detect 6.4 ±0.6 



multiple systems if all of the 39 stars in our sample were binaries. The successful detection of 9 multiple 
systems out of 39 AeBe stars corresponds to a mean number of companions of 1.4. Taking into account the 
uncertainty in our monte-carlo prediction and Poisson statistics we calculated that there is a 95% likelihood 



that the AeBe star binarity is greater than 85%. This is comparable to that of TTS ( 80% ±22% Leinert 



et al. (1993), phcz et al. (1993)| ) and exceeds that of near solar type MS stars (57%, pDuqucnnoy fc Mayor 



1991) ). Note that sources were previously detected around some of these AeBe stars (Z CMa ( Haas et 



al. (1993) ), VV Ser and V376 Cas ( Li et al. (1994) ) for example), but that these were not combined with 
our observations because they were not within the completeness region of this study (0.4"<sep.<8" and 
K <10.5). It should be noted that to obtain a significant increase in the confidence level would require a 
sample at least four times as large. Such a large sample is not available in the northern hemisphere, but 
this statistical study can be improved by observing additional AeBe stars in the southern sky. 

It is also important to note that we found companions (which are assumed to be the fainter stars) around 
30% of the type I AeBe stars in our sample (HD 200775, BD+46 3471, MWC 1080, KK Oph, and HD 
150193). In the case of HD 200775 and BD+46 3471 the companions were found to be significantly fainter 
( AK = 4.9 and 4.7 respectively) than the primary star and can only account for a small portion of the 
observed infrared light excess (priillcnbrand et al. (1994) ) . The nearer of the two companions around MWC 



1080 contributes about 10% of the 2.2 /im radiation and we estimate that the companions around KK Oph 
and HD 150193 contribute about 10% and 13% of the total 2.2 fim radiation. The companion near V892 
Tau, the only type II star with a companion in this study, is seen to contribute less than 1% of the total 2.2 



fira light. The three stars that have not previously been classified by Hillenbrand et al. (1994) have been 



found to all have very bright companions contributing between 20% and 100% of the total infrared light 
from the system. 

The ambiguity as to which star is the primary star and which one is the secondary cannot be resolved 
without further observations of these systems using different wavelengths. However, we can measure the K 
magnitude differences between the various components of the observed multiple systems and compare those 



to the 2.2 /im infrared excess estimates published by Hillenbrand et al. (1994). The measured AK of MWC 



1080, V892 Tau, and KK Oph are comparable to the estimates from Hillenbrand et al. (1994). The infrared 
excess in those systems could therefore be caused by the observed stars if the primary AeBe star is actually 
the dim star in those fields. The AK measured in the remaining three systems for which an estimate of the 
infrared excess is available are an order of magnitude too small and none of these companions are likely to 
be the source of the infrared excesses. No stars were found around other type I stars even in stars such as 
MWC 275 and AB Aur which are relatively close by and have been observed to have a large infrared excess 



(Hillenbrand et al. (1994)). We believe that the relatively close distance of some of these type I AeBe 
stars and the fact that no star was found close to them make them good candidates for the search for an 
extended structure. 



4. CONCLUSION 

We have successfully used the shift-and-add technique to obtain high resolution images of AeBe stars. 
We observed 39 AeBe stars and determined the binarity of AeBe stars to be high (85%) based on the 
detection of 9 multiple systems. While none of the companions seem bright enough to be a major source 
of infrared excess, we showed that the 2.2 /im infrared excess of three of the observed systems could be 
explained if we assume that the dim star is the primary AeBe star. The infrared excess in the remaining 
systems cannot be attributed to any of the observed companions. 

5. ACKNOWLEDGMENT 

Dr. E.J. Spillar would like to acknowledge the support of the US Air Force during his sabbatical at the 
Starfire Optical Range. We would like to thank Jeff Kuhn of Michigan State University and the National 
Solar Observatory and Steve Boughn of Haverford College for helping to build the camera. 



REFERENCES 

Abt, H. A., Levy, S. G., 1976, Ap. Suppl. 30, 273 

Allen, C. W., Astrophysical Quantities (London: Athlone) 

Bates, R. H. T., Cady, F. M. 1980, Opt. Commun., 32, 365 

Berrilli, F., Corciulo, G., Ingrosso, G., Lorenzetti, D., Nisini, B., Strafella, F. 1992, ApJ, 398, 254 

Blondel, P. F. C., Talavera, A., Djie H. R. E. T. A. 1993, A&A, 268, 624 

Brooke, T. Y., Tokunaga, A. T., Strom, S. E. 1993, AJ, 106, 656 

Christou, J. C. 1991, PASP, 103, 1040 

Duqucnnoy, A., Mayor, M. 1991, A&A, 248, 485 

Finkenzeller, U., Mundt, R. 1984, A&AS, 55, 109 

Fischer, D. A., Marcy, G. W. 1992, ApJ, 396, 178 

Ghez, A. M., Neugebauer, G., Matthews, K. 1993, ApJ, 106, 2005 

Haas, M, Christou, J. C., Zinnecker, H., Ridgway, S. T., Leinert, CH. 1993, A&A, 269, 282 

Hartmann, L., Kenyon, S. J., Calvet, N. 1993, ApJ, 407, 219 

Herbig, G. H. 1960, ApJS, 4, 337 

Herbig, G. H., Robbin-BeU, K. 1988, Lick Observatory Bulletin No. 1111 

Hillenbrand, L. A., Strom, S. E., Vrba, F. J., Keene, J. 1992, ApJ, 397, 613 

Hogbom, J. A. 1974, A&AS, 15, 417 

Kenyon, S. J., Whitney, B. A., Gomez, M., Hartmann, L. 1993, ApJ, 414, 773 

Leinert, Ch., Zinnecker, H., Weitzel, N., Christou, J., Ridgway, S. T., Jameson, R., Haas, M., Lenzen, R. 
1993, A&A, 278, 129 

Li, W., Evans, N. J. II, Harvey, P. M., Colome, C. 1994, ApJ, 433, 199 



Mamon, G. A., Soncira, R. M. 1982, ApJ, 255, 181 

Natta, A., Prusti, T., Krugel, E. 1993, A&A, 275, 527 

Terranegra, L., Chavarria-k, C, Diaz, S., Gonzalcz-Patino, D. 1994, A&AS, 104, 557 

The, P. S., de Winter, D., Perez, M. R. 1994, A&AS, 104, 315 

Wood, D. O. S., Myers, P. C, Daugherty, D. A. 1994, ApJS, 95, 457 



This manuscript was prepared with the AAS lAT^^X macros v4.0. 






H 


^ 


B 


P3 


Cp 


p 


X 


1 — 1 


fD 


^ 




&. 


y 


n 


O 


fD 


o 








o 


M 


o 


o 


d 


^ 


0) 


1 — ' 

02 








CO 




bO 




C3" 


CO 


O 


h^' 




t-t-j 


a 


M 


r1- 


o 


O 


Oi 


F? 


h^ 


I— ^ 


ri 


Ol 


1—^ 


n 


o 


a> 




fD 


1 — ' 

CO 


1—^ 

p 

r1 


P 


CO 








O 


f^ 






« 


'U^ 


(D 


O 


> 


X 




fD 

CO 

fD 


^Tl 


P3 

p 


CO 


fD 




►-i 


D- 


CO 


^. 


M 


fD 

B 


CO 


d 


fD 






3 

55' 






p 

P 


02 
CO 


to 

o 


p 




c_> 




H 

P" 

fD 




n 


O 


ra 


O 


a 


a 




o 


+ 


0! 


►-i 


02 


W 


, , 


CO 


rr 




>)^ 


?^ 


fD 


1 — ^ 


p 

p 


CO 


fef 


03 


s 


^ 


o, 


p 


O 


a- 


cr 


1 — ^ 


r+ 


r+ 


o 




►-i 


m 


c^ 


m 


o 




p 


< 

CX) 
CD 




a- 


to 



T31 

o 



€^ 




HD 200775 



arcsecond 



(cntrs=40,1 6,6,3,1) 




V892Tau (Ellas 1) 



-1 1 

arcsecond 

(cntrs=45,20,9,4,2) 



2 

§1 
O 

o 

-2 
-3 



•u 




KKOoh 



-3-2-10 1 2 3 

arcsecond 

{cntrs=1 00,40, 18,8,3) 



"D 



Cttl 



■ 'U 



BD-h46 3471 



-10 12 

arcsecond 

(cntrs=40,16,6,3,1) 



© 



XYPer 



-1 1 

arcsecond 

(cntrs=1 00,43, 18,7,3) 



b 
5 
4 

O o 


•o 'U 


-1 

-2 


m 


-3 

-4 


HD 141569 



-7 -6 -5 -4 -3-2-1012 

arcsecond 

(cntrs=40,26, 17,1 1,7) 



■D 

§^ 
O 

0)0 

o 



'U 




MWC 1 080a 




HD 216629 



'U 



© 



-2 -1 

arcsecond 

(cntrs=40,1 6,6,3,1) 



0.5 


'U 


c 
o 

o 

CO 


0^ 


•0.5 






MWC 1 080b 



-3-2-10123 

arcsecond 

(cntrs=32,10,3,1,0.3) 



-0.5 0.5 

arcsecond 

(cntrs=40,16,6,3,1) 



1 


^ 


C 

o 
oo 

U) 


jO 


CO 

-1 


^r 




HD 150193 



o 
i 



-1 1 

arcsecond 

(cntrs=100,35, 13,4,1) 



11 



Table 1. AeBe Stars 



Star 


Group^ 


D2 


Ref. 


q miii^ 


Disc. Rate* 






(pc) 


D 


(Af2/Mi) 


(%) 


VVSer 


I 


440 


6 


.2 


18 


MWC 614 








.2 


17 


WW Vul 




550 


7 


.35 


13 


V1295 Aql 








.15 


18 


AS 442 




700 


8 


.2 


17 


V1493 Cyg 








.35 


12 


HD 200775 


I 


600 


6 


.15 


19 


V645 Cyg 




6000 


12 


.2 


8 


V373 Ccp 


I 


1000 


5 


.3 


13 


BD+46 3471 


I 


900 


6 


.2 


16 


HD 216629 




725 


12 


.3 


14 


MWC 1080 


I 


2500 


6 


.3 


10 


V376 Cas 


H 


600 


6 


.3 


14 


VXCas 








.2 


17 


V594 Cas 


I 


650 


9 


.2 


17 


IP Per 




350 


10 


.4 


13 


XYPer 




160 


12 


.2 


22 


V892 Tau5 


II 


160 


5 


.2 


22 


AB Aur 


I 


160 


6 


.1 


26 


MWC 480 








.2 


17 


UXOri 


I 


460 


6 


.3 


15 


MWC 758 








.2 


17 


CQ Tau 








.2 


17 


HK Ori 


I 


460 


6 


.3 


15 


HD 244604 








.3 


14 


Z CMa 




1150 


9 


.1 


17 


MWC 166 


HI 


1150 


6 


.15 


16 


KK Oph 


I 


310 


12 


.2 


17 


HD 141569 








.3 


14 


HD 144432 








.37 


12 


HD 203024 








.35 


12 


LkHa 134 




700 


12 


.3 


14 


MWC 275 


I 


125 


11 


.2 


22 


MWC 297 


I 


450 


12 


.05 


23 



12 



Table 1 — Continued 



Star Group^ D^ Ref. q inin^ Disc. Rate* 

(pc) D (Ah/Ah) (%) 



MWC 300 




700 


8 


.2 


17 


HD 150193 


I 


160 


6 


.28 


14 


V1685 Cyg 


I 


1000 


8 


.1 


18 


V1686 Cyg 


I 


1000 


8 


.15 


16 


V361 Cyg 


III 


1000 


9 


.2 


15 



^As defined by Hillenbrand ot al. (1994) 



^Distance of object from the Hterature 

Minimum mass of a detectable companion. 
^Probability of detecting an existing secondary around the star. 



^5 V892 Tau listed as Elias 1 in [Hillenbrand et al. (1994) 
References. — (6) Hillenbrand et al. (1994) ;(7) Brooke et 



al. (1993)1 ; (8) [Terranegra et al. (1994)| (9) 



FinkenzcUer & Mundf 



(1984);(10) Wood et al. (1994) (11) Blondel et al. (1993);(12) 



BerriUi et al. (1992) 



13 



Table 2. Detected Multiple AeBe Systems 



Star 


Sep. 


Pos. Angle 


p.si 


AK2 


ARS 




(arcsec) 


(degrees) 


(A.U.) 






HD 200775 


2.25 


164 


1350 


4.9 


1.5 


BD+46 3471 


4.18 


43 


3762 


4.7 


1.0 


HD 216629 


6.96 


147 




0.0 




MWC 1080 


4.69 


86 


11725 


6.2 


2.76 




.6 


270 


1500 


3.1 


2.76 


XYPer 


1.2 


255 


192 


0.0 




V892 Tau2 


3.72 


19 


595 


4.2 


5.0 


KK Oph 


1.5 


257 


465 


2.5 


3.0 


HD 141569 


8 


314 




2.0 






6.8 


312 




2.7 




HD 150193 


1.1 


236 


176 


2.2 


5.0 



^Projected separation at distance D given in Tabic 1 

^K magnitude difference between primary and companion star 
^Froni Hillenbrand et al. (1994), infrared excess at 2.2 /im in 



astronomical magnitude 

*V892 Tau listed as Ellas 1 in Hillenbrand et al. (1994) 



